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ABSTRACT 

This  lecture  presents  the  general  guidelines  we  have  followed  to  build  up  an  exhaustive  and 
accurate  spectroscopic  database  for  radiative  transfer  in  air  and  CO2-N2  plasma  applications, 
including  atmospheric  entries.  This  High  Temperature  Gas  Radiation  (HTGR)  database  in¬ 
cludes  bound-bound  atomic,  diatomic  and  polyatomic  (in  particular  CO2)  transitions,  bound-free 
transitions  resulting  from  various  mechanisms,  and  free-free  transitions.  The  covered  spectral 
range  is  1000-200,000  cm -1  and  the  targeted  maximum  temperature  is  30,000  K.  A  particu¬ 
lar  attention  is  given  to  the  selection  of  the  most  reliable  fundamental  data  (line  strengths  for 
atoms  and  ions,  electronic  transition  moment  functions  and  spectroscopic  constants  for  diatomic 
molecules,  cross  sections  for  continuum  radiation),  and  to  accurate  prediction  of  atomic  and  ion 
line  shapes  which  are  strongly  self-absorbed.  The  spectroscopic  database  is  applied  in  the  last 
part  of  this  lecture  to  the  simulation  of  Fire  II  experiment  as  an  example  of  earth  entry,  and 
to  the  prediction  of  a  Mars  entry  simulation.  We  discuss  in  these  applications  the  spectral  dis¬ 
tribution  of  radiative  fluxes  from  the  infrared  to  the  VUV,  the  effects  of  optical  thicknesses  and 
self- absorption,  and  the  effects  of  chemical  nonequilibrium  on  continuum  radiation. 

1  Introduction 

During  atmospheric  entries,  the  gas  heated  in  the  shock  layer  will  emit  a  radiation  that  will 
contribute  to  the  incident  flux  at  the  vehicle  surface.  Different  mechanisms  are  involved  which 
give  rise  to  a  complex  structure  of  the  emission  and  absorption  spectra.  Bound-bound  transi¬ 
tions  in  atoms  and  molecules  will  give  line  spectra.  Free-free  transitions  will  contribute  to  a 
continuous  spectrum.  Bound-free  transitions  will  give  a  continuous  spectrum  with  eventually 
some  structures.  The  radiative  spectral  intensity  /CT(u)  (W.m-2.sr.(m-1)-1)  can  be  obtained 
at  each  point,  for  each  propagation  direction  u,  and  each  wavenumber  a  (nr-1)  by  solving  the 
radiative  transfer  equation  (RTE): 


dla(n) 

— - - =  rja  -  KaIa{u),  (1) 

as 

where  s  is  the  optical  path  along  u,  is  the  absorption  coefficient  (m-1),  and  r/CT  is  the  emission 
coefficient  (W.m-3.sr-1.(m-1)-1).  The  above  equation  has  been  written  assuming  that  scatter¬ 
ing  by  possible  solid  particles  is  negligible,  otherwise,  the  propagation  directions  are  coupled 
by  this  phenomenon.  The  RTE  requires  one  boundary  condition  corresponding  to  the  intensity 
leaving  the  wall  point  which  intersects  the  propagation  direction  if  u  is  directed  outward  this 
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wall,  and  zero  otherwise.  The  radiative  flux  on  a  surface  can  be  calculated  from  the  intensity 
following: 

da  /  I a 

J  47T 

where  n  is  the  normal  to  the  considered  surface,  and  the  amount  of  power  exchanged  between 
the  matter  and  the  radiation  field  per  unit  volume  is  given  locally  by: 


(u)  u  ■  n  du 


(2) 


Pr 


da  /  {r)0 


/  47 r 


KaIa{n))  du 


(3) 


Ka  and  r]a  are  local  quantities  which  depend  on  the  thernrochenrical  states  of  the  emitting  and 
absorbing  species.  If  the  medium  is  at  local  thermodynamic  equilibrium  (LTE)  at  temperature 
T,  Ka  and  r]a  are  linked  through  the  Kirchhoff’s  law: 

Va  =  KaI°a(T),  (4) 


where  J°(T)  is  the  equilibrium  intensity  given  by  the  Planck’s  law: 

2hc2a3 


m  = 


exp(  w)  —  i 


(5) 


h,  k  and  c  being  respectively  Planck  and  Boltzmann  constants,  and  the  speed  of  light. 


An  accurate  modelling  of  radiation  is  required  for  several  reasons: 

•  Radiation  is  currently  used  in  ground  experiments  and  in  some  flight  tests  to  characterize 
the  state  of  the  medium.  These  optical  measurements  are  usually  done  in  the  visible  and 
UV  range,  sometimes  in  the  IR  and  vacuum  UV  (VUV).  It  is  thus  necessary  to  have  data 
of  spectroscopic  quality  in  these  spectral  ranges. 

•  For  high  velocity  entries,  the  radiative  flux  at  the  wall  can  be  important  and  even  pre¬ 
dominant.  To  calculate  this  flux,  it  is  necessary  to  have  a  data  base  which  includes  all  the 
mechanisms  that  will  contribute  to  the  emission  or  absorption  in  a  wide  spectral  range. 

•  Radiation  is  an  additional  energy  transport  mechanism  and  can  also  affect  the  chemical 
state  of  the  gaz  through  emission,  photodissociation,  photionization  ,  radiative  recombina¬ 
tion  ...  In  some  situations,  it  can  be  necessary  to  take  into  account  the  effect  of  radiation  on 
the  flowfield  and  to  perform  coupled  calculations  of  the  fowfield  and  of  the  radiation  field. 
However  solving  the  RTE  is  computationally  very  expensive.  To  reduce  the  calculation 
cost,  approximate  models  of  radiative  properties  can  be  developed. 

There  are  several  computer  codes  and  spectral  databases  developed  for  radiation  analysis  (NEQAIR96 
[1],  LORAN  [2],  SPRADIAN  [3],  MONSTER  [4],  SPEC  AIR  [5],  PARADE  [6],  HARA  [7],  GPRD 
[8]. . . However  some  of  them  are  classified,  other  do  not  completely  satisfy  the  following  re¬ 
quirements: 

•  Cover  a  wide  temperature  range  (up  to  30000K  for  entry  applications  but  also  for  other 
plasma  applications). 
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Consequently,  cover  the  whole  useful  spectral  range  including  infrared  and  VUV  ranges. 

Allow  line-by-line  calculations  with  the  fine  structure  of  bound-bound  atomic  and  molec¬ 
ular  spectra  to  accurately  predict  radiative  transfer  with  absorption. 
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For  these  reasons,  we  have  developed  a  database,  called  HTGR  (High  Temperature  Gas  Radi¬ 
ation),  for  air  first  [9,  10]  using  the  latest  spectroscopic  constants  and  transition  probabilities. 
This  database  was  recently  extended  to  CO2GN2  plasmas  [11,  12]  and  the  development  of  some 
approximate  radiative  property  models  is  under  progress. 

In  the  following  we  give  the  guidelines  we  have  followed  to  build  this  database.  Such  devel¬ 
opments  rely  on  spectroscopic  concepts  which  are  not  straightforward.  The  interested  reader 
is  invited  to  consult  specialized  books  on  these  concepts.  A  few  examples  of  application  of  the 
database  are  given.  For  applications  involving  local  thermal  equilibrium,  we  calculate  the  pop¬ 
ulations  of  the  different  emitting  species  from  statistical  physics.  For  entry  applications  under 
nonequilibrium  conditions,  we  use  thermochemical  fields  calculated  by  other  groups. 

2  Database 

The  database  allows  to  calculate  the  monochromatic  emission  coefficient  rja  and  the  monochro¬ 
matic  absorption  coefficient  Ka  for  all  the  radiative  processes  and  species  listed  in  table  1  which 
are  relevant  to  Earth  and  Mars  entries.  Other  polyatomic  IR  data  (H2O  and  CH4)  have  also 
been  developed  for  combustion  applications  [13,  14]. 


Bound-bound  transitions 

-  atomic  species 

-  diatomic  species 

-  triatomic  species 

N,  N+,  N++,  0,  0+,  0++,  C,  C+ 

N2,  02,  NO,  N+,  C2,  CO,  CN,  CO+ 
C02 

Bound-free  transitions 

-  atomic  photoionization 

N  ,  N+,  0,  0+ 

-  atomic  photodetachment 

N“,  0",  C" 

-  molecular  photoionization 

N2,  02,  NO,  C2,  CO,  CN 

-  molecular  photodissociation 

O2  Schumann-Runge 

free-free  transitions 

N,  N+,  N++,  0,  0+,  0++,  C,  C+ 

other  mechanism 

chemiluminescence  CO  +  0 

Table  1:  Radiative  processes  taken  into  account. 


2.1  Bound-bound  transitions 

We  consider  a  radiative  transition  between  bound  energy  levels  u  (upper)  and  l  (lower)  of  a 
species  A.  Three  mechanisms  have  to  be  retained  for  energy  transfer  applications:  spontaneous 
emission,  absorption,  and  stimulated  emission,  corresponding  to  the  following  reaction  schemes 

A(u )  — >  A(l)  +  hccrui  spontaneous  emission 

A(l)  +  hccrui  — >  A(u)  absorption 

A(u)  +  hcaui  — >  A(l)  +  2hcaui  stimulated  emission 

where  crui  is  the  wavenumber  of  the  emitted/absorbed  photon.  hcaui  in  the  above  reactions 
designates  the  energy  of  the  photon  and  satisfies 

hcaui  =  Eu-  Ei,  (6) 

where  Eu  and  £)  are  the  energies  of  the  upper  and  lower  levels.  Three  Einstein  coefficients 
Auh  Biu  and  Bui  associated  respectively  to  the  three  above  radiative  mechanisms  are  usually 
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introduced,  enabling  to  express  the  contribution  of  these  bound-bound  transitions  to  the  spetral 
emission  and  absorption  coefficients  of  the  medium  at  current  wavenumber  a  according  to 

Vcr  ^  '  |  hcauinufui{cr  @u{), 

ul 

l^a  =  ^  ^uBul)h(Julful(&  °ul), 

ul 

where  nu  and  rq  designate  the  population  of  the  upper  ( u )  and  lower  ( l )  transition  levels,  and 
fui(cr  —  a ui)  is  the  spectral  line  shape  of  the  transition  accounting  for  Doppler,  natural  and 
collisional  broadening.  Einstein  coefficients  are  intrinsic  parameters  of  the  radiating  particle 
and  satisfy  the  relations 


(7) 

(8) 


Aui  =  8irhcahBui,  (9) 

9uBui  —  9i  Biu ,  (10) 

where  gu  and  gi  are  the  degeneracies  of  levels  u  and  l.  Because  of  these  relations  between  Einstein 
coefficients,  we  need  to  know  only  one  of  them  for  each  transition. 

To  predict  bound-bound  absorption  or  emission  spectra,  we  need  to  know  for  all  the  tran¬ 
sitions  the  line  position,  aui,  the  strength  of  the  transition,  e.g.  Biu  (in  fact  we  will  see  in  the 
following  that  atomic  and  molecular  spectroscopists  use  different  quantities  related  to  Biu),  the 
line  shape  fui(cr  —  crui),  and  of  course  the  thermodynamic  state  of  the  medium,  i.e.  the  nu  and 
ni  values. 

In  the  following,  we  will  discuss  only  electric  dipolar  transitions  since  they  are  the  main 
contributors  to  radiative  energy  transfer.  In  this  case,  the  absorption  Einstein  coefficient  Biu 
may  be  expressed  as 

B  —  1  hi 

J-Jlu  0J  9  .  -Oui 

3/Dc47T£0  gi 

where  the  transition  moment  Rui  is  obtained  from  the  the  dipolar  moment  operator  /i  of  the 
radiating  particle  according  to 


Rul  =  ^2  v\u,mu)\2  ■ 

mimu 


(12) 


The  summation  on  mi,  resp.  mu  in  the  above  expression  extends  over  a  complete  basis  of  the 


state  space  associated  to  the  level  l,  resp.  u. 
general  selection  rules: 

Pu 

A  J 
Ji  =  0 


Dipolar  radiative  transitions  satisfy  the  following 


+ 

Pi 

(13) 

= 

0,±1 

(14) 

Ju  =  0 

(15) 

where  pu  and  pi  are  the  parities  of  levels  u  and  l  (the  parity  of  a  state  designates  its  symmetry 
in  a  space  reflexion),  Ju  and  J;  are  their  associated  total  angular  momentum  quantum  numbers, 
and  A  J  =  Ju  —  Ji- 


2.1.1  Atomic  line  spectra 

A  one-electron  atomic  state  is  defined  by  the  quantum  numbers  n,  l,  mi,  s,  ms.  n  is  the  principal 
quantum  number  and  takes  positive  integer  value.  I  is  the  orbital  angular  momentum  quantum 
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number  and  takes  integer  values  between  0  and  n  —  1.  The  numerical  values  of  l  can  be  replaced 
by  letters  (s,p,  d,  f,g,h. . The  quantum  numbers  mi  and  ms  are  the  projections  of  the  orbital 
and  spin  momenta.  The  spin  quantum  number  s  is  equal  to  2 .  For  a  nrany-electron  atom, 
the  wavefunction  can  be  expressed  in  a  first  approximation  in  terms  of  product  of  one-electron 
atomic  states  nlmi  with  the  requirement  of  respecting  the  Pauli  exclusion  principle  (impossibility 
to  have  2  electrons  in  the  same  state,  i.e.  with  the  same  quantum  numbers).  This  allows  to  define 
a  configuration  by  counting  the  number  N  of  electrons  in  each  subshell  nl ,  using  the  notation 
nlN .  The  ground  state  configuration  for  atomic  nitrogen  is  denoted  by  ls22s22p5 ,  which  indicates 
that  two  electrons  are  in  the  subshell  Is,  two  are  in  the  2  s  one  and  three  are  in  the  2 p  one.  The 
two  first  excited  state  configurations  of  atomic  nitrogen  are  ls22s22p23s  and  ls22s2p4.  In  fact, 
we  have  to  take  into  account  the  coupling  of  the  different  angular  momenta  of  the  incomplete 
subshell  .  Two  extreme  cases  of  coupling  correspond  to  jj  and  LS  coupling.  The  jj  coupling 
occurs  when  the  leading  process  is  the  coupling  between  orbital  and  spin  angular  momenta  of  an 
electron.  In  the  LS  coupling  scheme,  which  corresponds  to  the  atoms  an  ions  considered  here, 
the  electron  orbital  angular  momenta  are  coupled  to  give  the  total  orbital  angular  momentum 
L.  the  electron  spin  angular  momenta  are  coupled  to  give  the  total  spin  angular  momentum 
S.  The  combination  of  a  particular  S  value  and  particular  L  value  gives  a  spectroscopic  term 
2S+1L  (the  numerical  values  of  L  =  0, 1, 2, . . .  are  replaced  by  letters  S,  P,  D, . . .).  The  total 
angular  momentum  J  is  the  sum  of  L  and  S.  The  spin-orbit  interaction  splits  the  energy  term 
2S+1L  into  2  x  min(L,  S)  +  1  levels  denoted  2S+1Lj  corresponding  to  the  various  total  angular 
momentum  values  J  =  \L  —  S\,\L  —  S\  +  1, . . . ,  L  +  S  allowed  for  the  term  2S+1L.  Each  energy 
level  2S+1Lj  has  a  2  J  +  1  degeneracy  (without  considering  the  nuclear  spin). 

The  main  radiative  bound-bound  transitions  observed  for  atomic  species  are  electric  dipolar 
ones  which  satisfy  the  general  selection  rules  given  by  Eqs.  (13-15). 

Moreover,  if  the  upper  and  lower  transition  levels  satisfy  the  pure  LS  coupling  scheme,  we 
have  the  additional  selection  rules  for  the  electric  dipolar  transition 


AS  =  0, 

(16) 

A  L  =  0,±1, 

(17) 

o 

II 

o 

II 

sq 

(18) 

where  we  have  introduced  the  notations  AS  and  A L  concerning  the  total  spin  and  total  orbital 
angular  momenta  of  the  electrons  in  a  similar  way  as  for  J. 

Different  spectroscopic  databases  are  available  which  give  for  each  transition  u  <->  l  the 
transition  frequency,  the  energy  of  the  upper  Eu  or  lower  level  Ei  of  the  transition,  the  level 
degeneracies  gu  and  (ji ,  and  the  Einstein  emission  coefficient  Au[  or  the  absorption  oscillator 
strength,  fiu  which  is  defined  according  to 


flu 


BluMul 

nre 

8tt2 mec  1  ,,,  ,  ,  >  |2 

-Vul—  2^  \{Lm\  ^\u,mu}\  , 


3  hq~ 


9l 


mimu 


(19) 

(20) 


where  re,  me  and  qe  are  respectively  the  classical  radius,  the  mass  and  the  charge  of  the  electron. 
Emission  oscillator  strengths  are  generally  defined  following  the  convention 

ful  =  --flu,  (21) 

9u 

and  even  more  symmetric  oscillator  strengths  often  are  given  in  databases  since  we  have  the 
relation 

9ifiu  9ufui ■  (22) 
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To  calculate  the  emission  and  absorption  coefficient  we  need  to  know  the  populations  nu  and 
ni .  At  LTE  the  population  of  level  i  is  given  by 


n  i  = 


giNa 

Q(T) 


exp 


(23) 


where  Q(T)  is  the  internal  partition  function  of  the  radiating  particle  at  temperature  T,  and  Na 
the  total  population  of  the  radiating  species.  For  a  multi-temperature  model,  we  use  Eq.  23  at 
the  temperature  of  electronic  excitation. 

For  air  plasma  we  have  compared  (cf  Ref.  [15])  three  available  databases:  (i)  the  NIST 
database  [16]  which  contains  critical  acclaimed  data  but  is  thus  not  a  priori;  (ii)  the  Hirata  and 
Horaguchi  [17]  database  which  has  been  obtained  by  merging  the  1989  updated  Kurucz’s  [18] 
one,  old  NIST  (1966,1988)  ones,  line  positions  from  the  Kelly  database  [19]  and  an  earlier 
(1987)  version  of  the  TOPABSE  data  base  (see  next  item);  (iii)  the  TOPBASE  database  [20,  21] 
which  is  purely  theoretical  and  does  not  include  the  fine  structure.  TOPBASE  is  an  atomic 
database  which  concerns  astrophysically  abundant  species  (atomic  number  Z=l-26)  and  has 
been  developed  within  the  international  Opacity  Project  since  1984  in  order  to  estimate  stellar 
envelope  opacities.  The  database  contains  the  most  complete  dataset  of  LS-coupling  term  ener¬ 
gies,  f- values  and  photoionization  cross  sections.  They  have  been  computed  in  the  close-coupling 
approximation  by  means  of  the  R-matrix  method  with  innovative  asymptotic  techniques  [22]. 
Two  requisites  were  claimed:  completeness  and  accuracy.  For  example,  the  database  includes 
1391  energy  LS  terms,  and  18037  electric  dipolar  transitions  concerning  the  neutral  oxygen  atom 
(01). 

Figure  1  shows  a  comparison  between  the  numbers  of  lines  included  in  the  above  cited 
databases  for  O,  N  species  below  150000  cm-1,  and  for  0+,  N"1-  species  below  300000  cm-1.  Note 
that  the  Hirata  database  includes  all  lines  from  Kelly  and  Kurucz  ones,  so  relevant  comparisons 
concern  HIRATA,  NIST  and  TOPBASE  databases.  It  must  be  again  emphasized  that  the 
TOPBASE  database  does  not  include  the  fine  structure,  so  the  numbers  of  lines  inlcuded  in 
TOPBASE  plotted  in  Fig.  1  are  in  fact  numbers  of  multiplets.  The  optically  thin  emission  of 
each  atom  or  ion  transition 


_  , 

^ul  —  A  hC(JuiTlu 
47T 


(24) 


has  been  calculated  from  these  three  databases  for  an  atmospheric  LTE  air  plasma  at  10000  K 
for  N  and  O  and  at  20000  K  for  N+  and  0+.  Figure  2  gives  an  example  for  N  and  a  tem¬ 
perature  of  10000  K.  Some  misplaced  multiplets  in  TOPBASE  are  observed  in  a  few  regions, 
which  normally  do  not  contribute  to  radiation.  A  few  discrepancies  exist  between  NIST  and 
TOPBASE,  with  transitions  lacking  in  NIST.  The  discrepancies  are  on  the  whole  larger  with 
Hirata  base.  The  NIST  database  has  been  found  to  be  complete  enough  to  predict  emission  for 
the  temperature  range  considered  here.  We  have  thus  retained  the  NIST  database  for  atom  and 
ion  line  transitions. 

Bound-bound  atomic  transitions  may  generally  not  be  assumed  to  be  optically  thin  in  many 
applications  (in  particular  for  resonance  transitions  involving  the  ground  energy  term  of  the 
atom).  We  thus  have  to  calculate  the  detailed  high  resolution  absorption  and  emission  spectra 
of  the  plasma  for  radiative  transfer  predictions.  The  spectral  line  shapes  fui(cr  —  crui)  introduced 
in  Eqs.  (7-8)  have  thus  to  be  determined.  An  extensive  litterature  has  been  devoted  to  the 
spectral  line  shapes  of  atomic  species;  some  recent  reviews  may  be  found  in  Refs.  [23,  24].  The 
main  mechanisms  which  have  to  be  accounted  for  in  the  case  of  atomic  species  are  Doppler 
and  collisional  broadening:  they  may  be  treated  independently  with  sufficient  accuracy  for  the 
prediction  of  radiative  transfer.  The  global  line  shape  will  thus  be  the  convolution  of  line  profile 
due  to  each  mechanism. 
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7000 


KELLY  KURUCZ  HI  RATA  NIST  TOPBASE 


Figure  1:  Number  of  O,  0+,  N  and  N+  lines  included  in  the  different  atomic  databases  below 
150000  cm^1  for  atomic  species  and  below  300000  cm'1  for  ions.  The  numbers  indicated  in 
black  correspond  to  atomic  lines  located  above  28500  cm-1  (from  Ref.  [15]). 


For  the  applications  considered  here,  radiating  particle  velocities  have  a  Maxwellian  distri¬ 
bution,  so  the  Doppler  profile  is  a  Gaussian 


fD{a  -  aul)  =  \l  —  —  exp 

7T  7d 


-In  2 


®  ®ul 
ID 


(25) 


whose  half  width  at  half  maximum  (HWHM)  7#  is  expressed  (in  the  same  units  as  aui)  as: 


7 D  =  Pul 


2  kTr  In  2 


mrcz 


(26) 


In  the  above  expression,  Tr  is  the  translation  temperature  of  the  radiator  and  mr  its  mass. 

Neutral  collisional  broadening  is  generally  accounted  for  under  the  impact  approximation 
(see  Ref.  [25]),  leading  to  a  Lorentzian  profile  of  HWHM  7/;, 


/l(c7  -  CFul) 


7 £1  1 

n  (7 ul)2+(a~aul)2' 


(27) 


7^  may  be  expressed  as  the  sum  of  a  van  der  Waals  contribution  7)[fw  (corresponding  to  an 
interaction  potential  proportional  to  d~6,  d  being  the  distance  between  particles)  and  of  a  res¬ 
onance  contribution  (potential  proportional  to  cG3)  from  the  usual  expressions  given  by 
Griem  [25]. 

The  van  der  Waals  contribution  is  obtained  from 


„,vdw 

lul 


V 


f  97rfi5|Ar2|  \  5 

^  16 m\El  J 


(28) 


where  Np  is  the  population  of  the  neutral  perturber  p,  Ep  the  energy  of  the  first  perturber 
excited  state  which  mainly  determines  its  polarizability.  A r2  =  r2  —  rf  is  the  difference  between 
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Figure  2:  Comparison  between  the  emissivities  of  O  lines  included  in  TOPBASE,  HIRATA  and 
NIST  databases  at  10000  K  (from  Ref.  [15]). 


the  mean  square  radii  of  the  radiating  particle  in  the  upper  and  lower  transition  levels  if  the 
perturber  and  the  radiator  are  particles  of  different  type;  for  identical  particles,  only  the  term 
r2  associated  to  the  level  u  or  l  of  the  same  parity  as  the  radiator  ground  level  is  included  in 
the  above  difference  —  the  other  level  will  contribute  to  the  resonance  broadening  (see  below). 
The  symbol  (•)  designates  an  average  according  to  the  distribution  of  the  relative  velocity  v  of 
the  perturber.  For  a  Maxwellian  distribution  at  temperature  Tp.  this  average  is  expressed  as 

'  2sT(f )  V  rnp  ) 

where  rnp  is  the  reduced  mass  of  the  radiator-perturber  couple. 

The  resonance  contribution  occurs  only  for  perturbers  of  the  same  type  as  the  radiator  and 
is  expressed  according  to 


(29) 


YuT 


167r2eomec2 


lju 


O', 


UJ 


+ 


(30) 


where  the  summation  extends  over  the  levels  j  of  the  radiator;  note  that  for  a  given  j  level  in 
this  summation,  only  one  of  the  two  terms  does  not  vanish  since  u  and  l  parities  differ. 

Last,  broadening  due  to  collisions  with  charged  particles  (electron  and  ions),  known  as  Stark 
broadening  has  to  be  accounted  for.  For  entry  plasmas  which  are  usually  at  quite  medium 
pressure  and  weakly  ionized,  and  for  non  hydrogenoid  species  considered  here,  the  impact  ap¬ 
proximation  may  be  used  for  electron  and  ion  contributions.  One  may  use  simple  expressions 
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obtained  from  semi-clasical  adiabatic  approach  [26]  which  are  similar  to  those  given  above  for 
neutral  contribution.  As  an  example,  for  a  neutral  radiator,  the  Stark  line  width  (HWHM)  will 
be  expressed  as  the  sum  of  a  dipolar  contribution 


^diP  =  — — 

lul  Vp2vrc 


(31) 


and  a  quadupolar  one 

'rT“  =  (32) 

where  the  AC 4  and  AC3  parameters  characterize  resp.  the  dipolar  polarization  interaction 
potential  (proportional  to  d~A)  and  quadrupolar  one  (proportional  to  d~ 3)  and  may  be  expreesed 
according  to  [26] 


ac4  = 

C.±u  —  C41, 

C4J  = 

i2h  ( 

f a0qP 

)2E^. 

/  e  aj£ 

7T mehc2  * 

K  Qe 

ac3  = 

\l  (ft?C 

i2M 

Birf)  -BuirlrfyJ 

2  Ih 
rne 


Qp 

<7e 


(33) 

(34) 

(35) 


In  the  above  expressions,  Ih  is  the  Rydberg  constant,  ao  the  Bohr  radius  and  qp  the  perturber 
charge.  Bj  and  Bui  constants  may  be  crudely  estimated  from  the  total  electronic  angular  mo¬ 
mentum  quantum  numbers  Lu  and  L/  of  the  radiator  according  to 


+  1  /  Lj 

15  V  0 
Bui  =  -2 BuBt  |  ^ 


2 

0 

U 


(36) 

(37) 


Similar  expressions  concerning  ionized  radiators  may  be  found  in  Ref.  [27]. 

More  sophisticated  semi-classical  approaches  [27,  28,  29,  30],  beyond  the  previous  adiabatic 
approach,  have  been  used  for  the  present  database.  They  enable  to  achieve  better  accuracy 
even  if  their  use  remains  quite  tedious  in  the  framework  of  radiative  transfer  applications  where 
systematic  calculations  are  required. 


2.1.2  Diatomic  line  spectra 

The  total  hamiltonian  Htot  of  an  isolated  diatomic  molecule  can  be  written,  without  taking  into 
account  the  hyperfine  structure  due  to  nuclear  spin  of  nuclei  [31],  as: 

2  h 2  Ne  h2 

H‘°‘ -  - E -E^v"  +  *4- + + y» + ft/  (^) 

0=1  “  i=i  e 

were  Ma  is  the  mass  of  nucleus  a  and  me  the  electron  mass.  The  first  two  terms  correspond 
to  the  nuclear  and  electronic  kinetic  energy.  The  following  terms  correspond  to  the  nuclear- 
nuclear  potential  energy,  the  nuclear-electron  potential  energy,  the  electron-electron  potential 
energy.  The  last  term  Hsf  includes  the  different  terms  linked  to  interaction  of  spin  and  angular 
momenta  leading  to  the  fine  structure: 

Hsf  =  HSo  +  HSs  +  HSr  (39) 
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Hso  =  A(r)L.S  is  the  so-called  spin  orbit  interaction  which  represents  the  interaction  of  the 
electron  spin  characterized  by  spin  angular  momentum  S  with  the  magnetic  field  created  by 
electron  motion  characterized  by  the  orbital  angular  momentum  L.  H$s  =  e(r)[35'?  —  S'2]  is  the 
so-called  spin-spin  interaction,  Sz  being  the  projection  of  S  on  the  internuclear  axis.  Hrr  = 
7(r)N.S  is  the  so-called  spin-rotation  interaction  where  N  is  the  sum  of  the  orbital  electronic 
angular  momentum  L  and  the  rotational  angular  momentum  R  (N=L+R). 

The  level  energies  and  wavefunctions  which  are  required  to  calculate  line  positions  (Eq.  6)  and 
transition  probabilities  (Eq.  12)  correspond  to  the  eigenvalues  and  eigenfunctions  of  the  total 
hamiltonian. 

Characterization  of  a  rotational  state 

The  Born-Oppenheimer  approximation  which  separates  the  electronic  and  nuclear  motions  al¬ 
lows  to  define  electronic  states.  Similarly  to  atoms  there  is  a  coupling  between  electronic  spin 
and  orbital  angular  momenta.  For  molecules  with  light  nuclei,  this  coupling  is  weak  and  we 
can  define  the  projection  of  the  total  electronic  orbital  momentum  L  on  the  intermolecular  axis. 
The  absolute  value  of  this  projection  is  noted  A.  An  electronic  state  of  a  diatomic  molecule  is 
characterized  by  the  label  2S,+1A  where  2 S  +  1  is  the  spin  multiplicity.  The  numerical  values 
0,1,2,. . .  of  A  are  replaced  by  X,  n,  A,  . . .  To  differentiate  states  with  the  same  S  and  A  values,  a 
letter  is  added  which  is  generally  X  for  the  ground  state,  then  A,B,C,. .  .for  a  spin  multiplicity, 

a, b,c,. .  .for  another  spin  multiplicity,  etc  ....  Additional  labels,  g  or  u,  can  be  added  to  char¬ 
acterize  for  homonuclear  molecules  the  symmetry  of  the  electronic  wave  function  with  respect 
to  the  inversion  operation  i  in  the  molecular  frame,  and  +  or  —  to  characterize  for  X  states 
the  symmetry  of  the  electronic  wavefunction  with  respect  to  the  symmetry  operator  av  which 
corresponds  to  a  reflexion  through  a  plane  containing  the  internuclear  axis. 

For  A  >  0  there  is  a  double  orbital  degeneracy  called  A-  doubling  which  will  be  lifted  in  the 
rotating  molecule.  The  electronic  degeneracy  is  thus  equal  to  (2  —  5o,a)(2 S  +  1). 

If  we  consider  the  rotating  and  vibrating  molecule,  a  rotational  level  (or  rovibronic  level)  will 
be  characterized  by  the  five  quantities  n,v,  J,i,p  which  correspond  to: 
n:  the  electronic  state 
v:  the  vibrational  number 

J:  the  total  angular  momentum  without  nuclear  spin  (J=  L  +  S  +  R) 
i:  the  spin  multiplet  component  which  takes  2S+1  values 

p:  the  parity  which  characterizes  the  symmetry  of  the  total  wavefunction  (without  nuclear  spin) 
with  respect  to  the  symmetry  operator  I.  The  parity  takes  (2-<5o,a)  values. 

The  multiplet  structure  will  depend  on  how  the  electronic  angular  momenta  couple  to  the 
molecular  rotation.  Several  limiting  coupling  cases  are  defined:  the  Hund’s  coupling  case  a, 

b,  c,  d,  . . .  We  are,  in  the  present  study,  concerned  by  the  cases  a  and  b.  In  Hund’s  case  a,  the 
interaction  between  the  molecular  rotation  and  the  electrons  motion  is  weak,  L  and  S  are  closely 
coupled  to  the  internuclear  axis,  their  projections  being  A  and  X.  Since  R  is  at  right  angle  to 
the  internuclear  axis,  the  projection  of  J  on  the  internuclear  axis,  fi,  corresponds  to  =  A  +  X. 
The  quantum  number  J  takes  the  values  |H|+1,  |H|+2,  ....  In  Hund’s  case  b,  S  does  not  couple 
to  the  internuclear  axis  and  recouples  to  N  =  L+R.  J  can  have  all  integer  values  from  N+S 
to  N-S.  An  intermediate  a/b  coupling  has  been  considered  for  most  of  the  states  considered  here. 


Energy  of  a  rotational  state 

The  energy  of  the  rotational  level  (n,v,  J,i,p)  can  be  expressed  as  the  sum  of  an  electronic 
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contribution,  a  vibrational  contribution  and  a  rotational  contribution: 

P )  =  Eei(ji)  +  Evib(n,  v )  +  Erot(n,  v,  J,  i,  p )  (40) 

A  very  crude  approximation  consists  in  considering  that  the  molecule  vibrates  like  an  harmonic 
oscillator,  rotates  like  a  rigid  rotor  and  to  neglect  the  fine  structure.  In  this  case  we  obtain  a 
very  simple  expression  for  the  energy: 

E^i{J,p)  =  Eel(n)  +  Hu;(v  +  l/2)  +  BJ(J  +  l)  (41) 

where  hui  is  the  vibration  frequency  and  B  is  the  rotation  constant. 


We  have  used  experimentally  determined  energies  or  more  precisely  experimentally  determined 
parameters  which  allow  to  recalculate  the  energies.  Theoretically  [32,  33],  the  complete  hamil- 
tonian  is  expressed  on  a  basis  set  of  rotational-vibration-electronic  wavefunctions  which  corre¬ 
spond  to  Hund’s  case  a.  A  series  of  contact  transformations  is  performed  to  remove  or  weaken 
off-diagonal  matrix  elements.  The  final  effective  or  model  hamiltonian  leads  to  energy  expres¬ 
sions  which  involve  a  limited  number  of  parameters  which  are  fitted  to  experimental  data.  We 
give  two  examples  of  such  developments,  one  for  a  simple  case,  the  singlet  state  1S,  and  one  for 
a  more  difficult  case,  the  doublet  state  2IIU.  The  expressions  used  for  the  other  types  of  states 
can  be  found  in  [15,  11]. 

The  level  energies  for  1S  state  are  given  by: 

El^(J)  =  Tv  +  Bvx  -  Dvx2  +  Hvx3  (42) 


with  x  =  J(J  +  1). 

The  reduced  matrix  for  the  2IIU  state  and  for  given  v ,  J,  and  parity  has  two  dimensions: 


M  = 


Ta  Ha 
Ha  T\ 


with  [34] : 


A 


TV)a (J)  —  Tv - — — b  ( Bv  —  Adv)x 2  —  Dv  [x4  +  x2  —  l]  + 


1  ±  x 


\Pv  +  Qv  (1  ±  x)\  +  HyX2  [; x 4  +  6.x2  -  7] 


TVil(J)  =  Tv  +  -E  +  {BV  +  ADv)  [x2  -  2]  -  Dv  [(x2  -  2)2  +  x2  -  1 
|  [x2  -  1]  +  Hvx2  [x4  +  7] 

HvA(J)  =  2 Dv  [x2  -  1]  2  -  [bv  +  V-E  +  ^  (1  ±  x)]  [x2  -  1]  3  + 

Hv  [x2  —  l] 2  [x2  (3x2  +  l)] 


(43) 


(44) 


(45) 


where  x  =  J  +  The  A-doubling  is  accounted  for  through  the  ±  values  in  these  expressions. 


The  diagonalization  of  this  matrix  leads  to  the  following  expressions  for  the  energies  for  the 
two  multiplet  components: 
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Radiative  transitions 

A  rotational  line  will  correspond  to  the  transition  between  two  rotational  states  characterized  by 
(n7,  v' ,  J' ,  i',p')  and  (n7/,  v" ,  J" ,  i" ,p")  such  that  the  matrix  element  of  the  electric  dipole  moment 
operator  (Eq.  12)  does  not  vanish.  Such  transitions  obey  to  the  general  selection  rules  given  by 
Eqs.  (13-15)  and  to  additional  selection  rules  for  Hund’s  cases  a  and  b: 

AA  =  0,±1,  (48) 

AS  =  0.  (49) 

The  set  of  transitions  between  rovibrational  levels  of  the  electronic  states  A  and  B  will  be  called 
the  A-B  system.  The  transitions  between  the  vibrational  level  v’  of  A  and  the  vibrational  level 
v”  of  B  will  be  called  the  v’-v”  vibrational  band.  Finally  transitions  between  A,v’  and  B,v”  with 
the  same  AN  and  A  J  are  called  branchs  and  are  represented  by: 

YXlj  (50) 

where  i  and  j  designate  the  substates  of  the  spin-multiplet  of  the  upper  level  and  of  the  lower 
level.  X  and  Y  correspond  to  the  change  in  J  and  N  and  are  denoted  by  the  symbols  given  in 
Fig.  3.  The  different  branchs  corresponding  to  a  2IIr  — 2  nr  system  are  shown  schematically  in 
Fig.  4. 
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Figure  3:  Notation  of  the  different  branches 
according  A  J  and  AN  values. 


N+2 


M=-2  AW=-1  AN=0  jiN=l  jiN=2 


Figure  4:  Branch  strucure  for  a  2IIr  — 2  IIr 
system  when  the  two  states  are  close  to 
Hund’s  case  b;  N  is  even  in  this  example 
(from  [9]). 


The  strength  of  the  different  transitions  will  depend  on  the  transition  probability  Ru[  defined 
in  Eq.  12  which  is  written  as  the  product  of  two  terms 


Rui 


( Kv  fs. 


V  V 

J'J"' 


(51) 


(i?g,^',,)2  is  the  square  of  the  electronic  vibrational  transition  moment,  and  Svj,vj„  is  the  Honl- 
London  factor.  These  quantities  are  defined  following  the  convention  of  Whiting  et  al.  [35]  which 
results  in  the  uniform  sum  rule  for  the  Honl-London  factors: 


]  <S j'Vj"  —  (2  —  <5o,a,^o,A")(2S'  +  1)(2 J  +  1)  (52) 

r 
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The  electronic-vibrational  part  of  line  intensity  is  given  by: 


2 


il)v>(r)Re(r)i/)vii(r)dr  , 


(53) 


where  r  is  the  internuclear  distance,  Vv(r)  and  VV'(r)  are  the  radial  rotationless  vibrational 
wavefunctions  in  the  vibrational  levels  v'  (of  the  upper  electronic  state  of  the  transition)  and  v" 
(of  the  lower  electronic  state)  respectively.  I?e(r)  is  the  electronic  transition  moment  function 
(ETMF). 

We  need  to  know  for  all  the  contributing  bands.  This  information  can  be  extracted 

from  spectroscopic  measurements  but  is  far  to  be  complete.  To  remedy  this  lack  of  experimental 
data,  the  electronic-vibrational  part  of  line  intensities  have  been  calculated  systematically.  The 
method  comprises  three  steps. 

In  a  first  step,  the  potential  curves  for  all  the  elecronic  states  have  been  reconstructed  via  the 
Rydberg-Klein-Rees  (RKR)  procedure  [36,  37,  38]  which  allows  to  calculate  the  classical  turning 
points  r inner (v),  router{v )  from  the  vibrational  term  energy  G(v)  and  rotational  spectroscopic 
constant  Bv  through  the  equations: 


where: 


Tinneriv)  -  \j  +  f(V)  f(V )) 

(54) 
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In  Eqs.  (56,57),  fir  is  the  reduced  mass  of  the  molecule,  vq  is  the  zero  of  G(v),  G(v)  and  Bv  are 
Dunham  expansions  in  (v  +  ]>): 

n 

G(v)  =  YlYi0 
i= 0 
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The  resulting  potential  curve  Urkr{t )  was  extrapolated  by  using  a  repulsive  function  Uinner(r ) 
(inner  turning  points)  and  a  Hulburt  and  Hirschfelder  potential  [39]  Uouter(r )  (outer  turning 
points): 
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where  De  is  the  dissociation  energy  of  the  electronic  state  considered  relative  to  the  bottom  of 
the  potential  energy  curve,  and  re  is  the  equilibrium  internuclear  distance,  ai,  b\  and  02,  &2; 
C2  and  re  are  the  floating  parameters  adjusted  in  order  that  each  extrapolation  function  fits  a 
significant  part  of  each  extremity  of  the  RKR  potential  Urkr(v). 

In  a  second  step,  the  resulting  potential  Ufinai(r)  is  used  in  the  radial  Schrodinger  equation  to 
solve  for  the  rotationless- vibrational  wavefunctions  i/)v(r)  using  a  Chebyshev  expansion  method. 
We  checked  that  the  difference  between  reconstructed  and  initial  energies  remains  small  (lower 
than  0.14%  on  the  maximal  vmax). 

Finally,  (^R^,v"^j  is  evaluated  by  using  electronic  transition  moment  function  Re(r )  selected 
in  the  litterature  [9,  12].  Figure  3  shows  an  illustration  of  the  above  described  procedure  in  the 


Figure  5:  Reconstructed  potentials  for  the  O2  Schumann- Runge  (B-X)  system.  The  RKR  re¬ 
constructed  potentials  are  represented  by  points.  A  few  vibrational  functions  are  represented 
with  the  selected  transition  moment  function  (adapted  from  [9]). 

case  of  the  O2  Schumann-Runge  -  X3S~  transition. 

Diatomic  database 

We  have  calculated  an  exhaustive  spectroscopic  database  including  the  main  electronic  systems 
of  the  diatomic  molecules  N2,  N^~,  O2,  NO,  CO,  C2,  CN,  CO+  which  contribute  to  radiative 
transfer  in  N2  -  O2  -  CO2  plasmas.  These  systems  are  listed  in  Table  2.  For  each  of  these 
35  systems,  potential  curves  have  been  calculated  according  to  the  above  described  procedure 
from  up-to-date  deperturbed  —  when  available  —  Dunham  coefficients.  Electronic  transition 
moments  have  been  critically  selected  in  the  literature  and  squared  electronic  vibrational  tran¬ 
sition  moments  have  been  systematically  obtained  for  the  bands  listed  in  Table  2.  Systematic 
validations  have  been  achieved  by  comparing  experimental  or  theoretical  vibrational  radiative 
life  times  available  in  the  literature  to  those  calculated  from  our  vibrational  transition  moments. 
Rotational  line  positions  and  strengths  have  thus  been  calculated  from  up-to-date  carefully  se¬ 
lected  spectroscopic  constants.  Fine  structure  has  been  introduced  considering  the  intermediate 
Hund’s  case  between  a  and  b  cases  (E  —  E  transitions  have  been  calculated  considering  the  pure 
Hund’s  case  b).  A-doubling  has  been  accounted  for  when  data  were  available.  Honl-London 
factors  of  Ref.  [40]  have  been  used.  Regarding  perturbations,  rotational  perturbations  which 


8  - 14 


RTO-EN-AVT-162 


Radiation  Database  for  Earth  and  Mars  Entry 


Table  2:  Calculated  diatomic  electronic  systems 


molecule 

system 

upper  state  - 

Calculated  bands 

name 

lower  state 

(0  ■  0  ^ max ) 

n2 

First-Positive 

B3ug-A3E+ 

(0  :  21;  0  :  16) 

Second-Positive 

C3Uu~B3Ug 

(0  :  4;  0  :  21) 

Birge-Hopfield  1 

61nn-Y1S+ 

(0  :  19;  0  :  15) 

Birge-Hopfield  2 

6'1S+-Y1S+ 

(0  :  28;  0  :  15) 

Carroll- Yoshino 

c'41S+-V1S+ 

(0  :  8;  0  :  15) 

Worley- Jenkins 

c31ilu-x1x+ 

(0  :  4;  0  :  15) 

Worley 

031nu-v1s+ 

(0  :  4;0  :  15) 

NJ 

Meinel 

A2nu-x2z+ 

(0  :  27;  0  :  21) 

First-Negative 

b2e+-x2e+ 

(0  :  8;0  :  21) 

Second-Negative 

(0  :  6;  0  :  21) 

NO 

7 

A2E+-x‘2nr 

(0  :  8;  0  :  22) 

P 

B2  Ur-X2Ur 

(0  :  37;  0  :  22) 

5 

c2ur-x2ur 

(0  :  9;  0  :  22) 

e 

D2T,+-X2Ur 

(0  :  5;  0  :  22) 

i 

E2Z+-X2  Ur 

(0  :  4;0  :  22) 

P' 

B,2A-X2Ur 

(0  :  6;  0  :  22) 

11000  A 

d2e+~a2e+ 

(0  :  5;  0  :  8) 

infrared 

X2Ur-X2Ur 

(0  :  22;  0  :  22) 

o2 

Schumann-Runge 

B3E--X3E- 

(0  :  19;  0  :  21) 

CO 

Infrared 

Y1S+  -  X1^ 

(0:50;0:50) 

Fourth  Positive 

AlB  —  X3E+ 

(0:23;0:50) 

Hopfield-Birge 

R1S+  —  Y1S+ 

(0:2;0:50) 

Third  Positive 

63S+  -  a3n 

(0:2;0: 18) 

CO+ 

Comet-tail 

A2Ui  -  X2S+ 

(0:33;0:31) 

Baldet- Johnson 

r2s+  -  A2Hi 

(0:33;0:50) 

First  Negative 

R2S+  -  V2S+ 

(0:22;0:35) 

CN 

Red 

A2  Ui  -  V2E+ 

(0:38;0:34) 

Violet 

R2S+  -  V2S+ 

(0:25;0:36) 

LeBlanc 

r2s+  -  a2it 

(0:25;0:38) 

c2 

Philips 

A2Uu  -  WS+ 

(0:35;0:21) 

Mulliken 

£>1E+-X1E+ 

(0:22;0:21) 

Deslandres-D’Azambuja 

ClAg  -  A1  Uu 

(0:9;0:32) 

Ballik  et  Ramsay 

b3T,~  -  a3  PR 

(0:41;0:39) 

Swan 

d3n9  —  a3Iiu 

(0:18;0:33) 

Fox-Herzberg 

e3n g  -  a3nu 

(0:15;0:35) 
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affect  a  few  J  values  of  a  given  vibrational  level  have  not  been  modelled.  For  strong  homoge¬ 
neous  perturbations  such  as  those  encountered  for  the  five  N2  bnYj+,  03*^,  and 

excited  states,  and  for  the  two  NO  B2 II,.  and  C2Iir  excited  states,  perturbed  molecular 
constants  have  been  used. 

Rotational  lines  have  been  calculated  up  to  the  last  rotational  quantum  number  below  the 
centrifugal  dissociation  of  the  upper  and  lower  vibrational  levels,  i.e.  widely  beyond  the  validity 
range  of  the  selected  spectroscopic  constants.  Such  extrapolations  are  required  insofar  as  this 
database  is  devoted  to  predict  radiative  transfer  inside  very  high  rotational  temperature  plasmas. 
We  have  ensured  that  no  aberrant  line  positions  occur  due  to  these  extrapolations. 

References  from  which  the  various  selected  spectroscopic  constants  and  electronic  transition 
moments  have  been  taken  are  given  in  Ref.  [9]  for  N,  O  plasma  molecules  and  from  Ref.  [12]  for 
carbonaceous  molecules.  Detailed  systematic  validations  for  each  system  may  also  be  found  in 
these  two  references. 

For  each  of  these  lines,  a  Voigt  line  profile,  taking  into  account  Doppler,  collisional  and, 
when  required,  predissociative  broadening  has  been  adopted.  For  IR  rovibrational  transitions, 
detailed  descriptions  of  collisionnal  broadening  parameters  are  available  in  Refs.  [41,  42]  for  NO 
and  in  Ref.  [43]  for  CO  IR  systems.  However,  data  concerning  collisional  broadening  of  diatomic 
rovibronic  lines  remain  quite  scarce.  We  have  used  for  all  NO  UV  lines  the  collisional  HWHM 
from  Ref.  [44] 

jNOuv  —  0.293  f  J  cm_1atm_1,  (62) 

where  T  is  the  translational  temperature  of  the  heavy  particles.  For  all  other  systems,  we  have 
used  the  pragmatic  collisional  line  width 


7c  =  0.1 


(63) 


following  the  values  obtained  in  Ref.  [45] .  Finally,  we  have  introduced  predissociative  (lorentzian) 
line  widths  form  Refs.  [46,  47,  48,  49,  50,  51],  [52],  and  [53]  for  the  N2  VUV,  NO,  and  O2  systems, 
respectively. 

The  database  enables  to  calculate  any  absorption  or  emission  spectra  for  the  considered 
electronic  systems  from  the  populations  of  molecular  energy  levels  according  to  Eqs.  7-8.  For 
media  at  LTE,  the  populations  are  caculated  in  a  straightforward  manner  from  the  Boltzmann 
distribution  and  from  the  chemical  composition.  For  nonequilibrium  media,  a  currently  used 
approximation  is  the  two-temperature  model,  Trot  for  the  translational  motion  of  heavy  particles 
and  rotation  of  molecules,  Tve  for  molecular  vibration,  electronic  excitation,  and  free  electron 
translation.  In  the  framework  of  this  two-temperature  model,  the  population  of  an  energy  level 
(n,  v,  J,  i,  p )  is  expressed  as 


to. 


nvJip 


=  n7 


9snJip 


2  j  + 1 


Q(TVI 


,  Trot ) 


exp 


Eei(n)  +  Evib(n ,  v)  Erot(n,  v,  J,  i,p) 


kTn 


kTr 


rot 


(64) 


where  Q(Tve.  Trot)  is  the  two-temperature  partition  function  [54],  nr  is  the  total  population  of 
the  considered  molecule.  The  factor  gSnjip  is  the  nuclear  spin  statistical  factor  which  applies 
only  for  homonuclear  molecules. 


2.1.3  CO2  spectra 
Infrared  contribution 

It  is  well  known  that  CO2  is  very  active  in  the  infrared.  In  its  fundamental  electronic  state, 
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the  molecule  is  linear  and  possesses  three  modes  of  vibration:  a  symmetric  stretching 
mode  v i,  a  bending  mode  V2  (doubly  degenerated)  and  an  asymmetric  stretching  mode  1/3.  The 
three  modes  are  close  to  1351,  672  and  2396  cm-1.  The  modelling  of  CO2  radiative  properties 
at  high  temperature  is  a  challenging  problem  due  to  the  different  couplings  which  have  to  be 
correctly  described:  Fermi  resonance  between  symmetric  and  bending  modes  which  gathers  the 
levels  in  polyads,  Coriolis  effect  which  couples  bending  and  antisymmetric  modes,  ... 

The  spectroscopy  of  CO2  in  the  infrared  has  been  the  subject  of  a  lot  of  experimental  studies. 
The  HITRAN  [55]  molecular  spectroscopic  database  compiles  the  up-to-date  spectroscopic  data 
for  this  molecule  and  38  other  molecules.  It  is  regularly  improved.  However  this  database 
has  been  primarily  developed  for  low  temperature  applications  and  is  not  complete  at  high 
temperature;  it  does  not  include  the  transitions  from  high-excited  rovibration  states. 

A  few  attempts  have  been  made  to  generate  data  suitable  for  use  at  high  temperature  especially 
for  gases  involved  in  combustion  and  in  particular  for  CO2.  The  HITEMP  [56]  database  has  been 
calculated  using  a  direct  diagonalization  method  and  is  adapted  to  temperatures  up  to  1000  K.  A 
similar  technique  but  with  a  better  treatment  of  the  Coriolis  operator  has  been  used  to  generate 
a  basis  for  applications  up  to  2000  K  [57] .  An  approach  based  on  effective  operators  has  been 
used  by  Tashkun  et  al  to  generate  the  CDSD-1000  database  [58]  which  is  able  to  reproduce 
successfully  the  high-resolution  spectra  at  low  and  medium  temperatures  and  satisfactorily  the 
medium  resolution  spectra  at  temperatures  up  to  3000K.  We  use  the  HITELOR  database  [59]. 
This  database  enlarges  the  HITRAN  data  by  using  approximate  and  empirical  formulas  for 
vibrational  band  intensities  and  by  computing  energies  through  a  diagonalization  of  the  Fermi 
matrix  associated  to  the  different  polyads.  The  HITELOR  databank  has  been  succesfully  used 
to  simulate  medium-resolution  spectra  at  temperatures  up  to  2850  K.  The  database  covers  the 
2.7,  4.4  and  15  jim  regions  which  are  the  most  important  in  terms  of  radiative  transfer.  All  the 
transitions  associated  to  a  lower  vibrational  level  lower  than  13000  cm-1  have  been  accounted 
for.  The  maximal  rotational  number  J  is  200. 

The  pressure-broadening  parameters  and  their  temperature  dependance  have  been  calculated 
semiclassically  up  to  J=101  [60].  For  J  >  101,  the  broadening  coefficients  have  been  set  to  their 
J=101’s  values. 

Figures  6,  7,  8  and  9  show  the  calculated  absorption  coefficients  for  a  C0(50%)-C02(50%) 
mixture  at  3600  K  and  a  pressure  of  100  Pa.  These  values  are  typical  of  those  encountered  in  a 
Mars  entry  (see  Sec.  3). 


VUV,  UV  and  visible  contributions 

As  could  be  expected  the  electronic  spectra  of  CO2  is  rather  complicated.  Two  phenomena 
have  been  retained:  the  UV  and  VUV  absorption  and  the  chemiluminescence. They  have  been 
observed  in  shock  tubes  and  flames  and  could  have  a  potentially  significant  contribution  in  the 
considered  applications. 

The  absorption  cross  sections  of  CO2  at  room  temperature  are  available  in  the  spectral  range 
50000  -  83000  cm”1  [61,  62].  They  have  also  been  measured  in  shock  tubes  between  28600 
and  52000  cm”1  up  to  3050  K  [63,  64],  In  this  spectral  region,  the  absorption  is  spectrally 
smooth  and  is  mainly  attributed  to  predissociative  transitions  in  the  X1Yj+  — A  Bo  system.  The 
contribution  of  this  mechanism  has  been  estimated  by  using  the  analytical  expression  of  the 
cross  section  o-(A,  T)  proposed  in  Ref.  [64]  : 

°(\T)  =  7p^exp(-ei(A)/UT)  (65) 

Wvib 

where  Qvib  is  the  vibrational  partion  function  of  CO2  and  ei(A)  corresponds  to  an  effective 
absorbing  level. 
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CO  :  Av=l 

100  Pa,  3600 K 
50%  CO 


!  0.2 


1600  1800  2000  2200  2400 

wavenumber  (cm'’) 

Figure  6:  Absorption  coefficient  for  the 

Av=l  bands  of  CO. 
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Figure  8:  CO2  absorption  coefficient  in  the 
4.3  /xrn  region 
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Figure  7:  CO2  absorption  coefficient  in  the 
15  fii n  region 

0.005 


CO, :  0«”o  ->  l«"l  <2.7  unit 
0.004  1  ' 


wavenumber  (cm*') 

Figure  9:  CO2  absorption  coefficient  in  the 
2.7  fim  region 


The  chemiluminescence  resulting  from  the  radiative  combination  of  carbon  monoxide  and  atomic 
oxygen, 

CO  +  O  — >  CO2  +  hca,  (66) 

is  also  believed  to  involve  the  electronic  transition  1B2~  >  X 1S+.  The  emission  coefficient  may 
be  expressed  as: 

Vo-  =  NoNcoAv  exp (ba/T)  (67) 

where  No  and  Noo  are  the  populations  of  O  and  CO.  The  parameters  Aa  and  b, T  have  been 
adjusted  on  the  shock  tube  data  of  Slack  and  Grillo  [65,  66]  in  the  spectral  range  15000  -  40000 
cm-1  and  for  temperatures  ranging  from  1300  to  3000  K  [11]. 

2.2  Bound-free  transitions 

Bound- free  transitions  correspond  to  photonization,  photodissociation  and  photodetachment. 
They  can  be  represented  formally  by: 

AB{k )  +  hca  — ►  A(i)  +  B(j)  (68) 
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The  energy  conservation  leads  to: 

hca  =  ^gg2  +  Ef  +  Ef  -  EfB,  (69) 

where  /i  is  the  reduced  mass  of  the  the  products  (A,B),  g  is  the  relative  velocity  of  the  products, 
Ef  is  the  energy  of  level  l  of  species  X.  The  emission  coefficient,  the  induced  emission  coefficient 
and  the  true  absorption  coefficient  can  be  expressed  as  a  function  of  the  level-to-level  cross 
sections  relative  to  the  spontaneous  emission,  Sf?k(a),  to  the  induced  emission,  S?j k(a),  and  to 
the  absorption,  Sfffa).  by: 


Vtrjb 

=  n?nf  hca  9sijAa)47T92f0{g)  ^ 

ijk 

(70) 

Kei 

=  hcagSflk{a)ATrg2f°(g)(^^- 

ijk 

(71) 

ahs 

Ka,fb 

=  Y.niBhmSi%(A-r 

(72) 

ijk 


where  nf  designates  the  population  of  the  level  l  of  species  X,  and  f°(g)  is  the  distribution 
function  of  the  relative  velocities  of  the  dissociation  products. 

By  considering  that  at  equilibrium  the  radiation  intensity  must  follow  the  Planck’s  law,  we 
can  write  relations  similar  to  the  Einstein-Milne  ones  [67]  between  the  three  level-to-level  cross 
sections: 

2 hca3  (73) 

(74) 

Assuming  a  Maxwell-Boltzmann  distribution  at  Trei  for  the  relative  velocity  distribution  function 
enables  to  write  the  emission  and  absorption  coefficients  as: 


qes 

Jij,k 


O) 


cei 

Dij,k 


{<T 


2  „2 


4t rg2g 
K 2 


9i9jStlk(°)  = 


Vajb 


2/mV  ^2 


A  R 

nfnf 


9k 


h°St-Aa)  exp 


ijk 


Trei)  gigj 


'Ef  +  Ef  -  EfB  -  hca ' 


kTr 


rel 


ijk 


h°St%{a)  nk 


AB 


1  - 


9k 


A  R 

nfnf 


Ef  +  Eb  -  Ef*  -  hca 


<AB 


9i9j  £(/b  Trei)n 


AB  exP 


kTr 


rel 


(75) 

,(76) 


where  g%  is  the  degeneracy  of  level  i  and  £(g,Trei)  is  the  volumetric  translational  partition 
function  of  a  particle  of  mass  g  at  temperature  Trei 


(  2irgkTrei  \  2 

V  h 2  ) 

Similarly  to  bound-bound  transitions,  the  exact  treatment  of  bound-free  transitions  in  nonequi¬ 
librium  requires  to  know  detailed  level-to-level  cross  sections,  which  are  not  always  available. 
We  detail  in  the  following  the  status  for  the  different  bound-free  mecanisms. 
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2.2.1  Atomic  photoionization 

The  atomic  photoionisation  corresponds  to  : 

A(k)  +  hca  — ►  A+(i )  +  e  (78) 


The  TOPBASE  database  [20]  provides  the  absorption  cross  sections  S&  from  levels  up  to  a 
principal  quantum  number  equal  to  10  or  11  as  a  function  of  photon  wavenumber.  These  cross 
sections  correspond  to: 

6k(a)  =  haY^Stfij  (79) 

ij 

where  j  corresponds  to  the  double  degenerate  unique  internal  energy  level  of  the  electron.  In 
the  two-temperature  model,  Eqs.  75,  76  and  79  lead  to  express  the  emission  and  absorption 
coefficients  as  : 


Va 


xneq 


na  t 

QaX 


Y  9k  exp 


L  k 

Ikion'IT-e— 


kTve 


<3fc(c) 


1  -xneqexp 


Q  at 


flat  2Qiont;(me,  Tve) 


exp 


Eu 


kTv, 


(80) 

(81) 

(82) 


where  nat,  n\ on,  and  ne_  are  the  total  populations  of  respectively  the  atom,  the  produced  ion 
and  electron,  Qat  and  Qmn  are  the  atom  and  ion  partition  functions  and  me  is  the  electron  mass. 
Eion  is  the  ionization  energy  corrected  from  the  Debye  ionization  lowering. 

The  factor  xneq  characterizes  the  deviation  from  equilibrium  and  consequently  from  Kirchhoff’s 
law. 


2.2.2  Molecular  photodissociation 

We  have  included  in  the  database  the  Schumann-Runge  continuum: 

02{k)  +  hca  — >•  0(i)  +  O(j)  (83) 

The  thermal  absorption  cross  sections  for  Schumann-Runge  continuum  are  available  in  the  tem¬ 
perature  range  300-10000K  and  the  wavenumber  range  48000-77000  cm-1  [68,  69].  Since  to  our 
knowledge  no  detailed  level-to-level  absorption  cross  sections  are  available,  we  need  to  make 
some  assumptions  to  deduce  the  absorption  and  emission  coefficient  in  the  framework  of  the 
two-temperatures  model  [70].  If  we  assume  that  the  oxygen  atoms  are  only  produced  in  the  3P 
and  1D  term  of  the  oxygen  ground  configuration  which  corresponds  to  the  asymptote  state  of 
the  O2  electronic  state,  the  thermal  cross  section  can  be  written  as: 

eLTE(a:  T]  =  qBt]  ^  9k  exp  (“  fr  )  (84) 

where  iq  and  jo  represent  the  3P  and  lD  terms  of  atomic  oxygen  and  Qo2(T)  is  the  internal 
partition  function  of  O2.  This  equation  combined  to  Eqs.  75  and  76  allows  us  to  deduced  the 
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spontaneous  and  induced  emission  coefficients  in  the  two-temperatures  model: 

Vcr  =  2/icV3exp  6LTE(a,Trot)n02Xneq 

C  =  Xneqno2 SLTE(ct,  Trot)  exp 

neq  _  rf) _ Qo2(Trot)  g  ( -^diss  +  E°D  _  E°D  \ 

no2Q20(Tve)t>(mo,Trot )  ^  kTrot  kTve  J 

where  no  and  no2  are  the  atomic  and  molecular  oxygen  total  populations,  -Ediss  the  dissociation 
energy  of  the  ground  electronic  state  of  O2,  and  E®D  the  energy  of  the  lD  term  of  atomic  oxygen. 
The  true  absorption  coefficient  is  approximated  by: 

nfs  =  n026LTE(a,Tve)  (88) 

%neq  in  the  above  equations  deviates  from  unity  in  chemical  nonequilibrium  which  leads  to  a 
deviation  from  Kirchhoff’s  law. 


(85) 

(86) 

(87) 


2.2.3  Molecular  photoionization 

Most  of  the  available  data  for  this  mechanism  correspond  to  room  temperature.  As  its  contribu¬ 
tion  remains  small  in  the  considered  applications,  we  use  pragmatically  the  room  temperature 
cross  sections  E  for  O2,  N2,  CO,  C2  and  CN  [71,  72]  and  evaluate  the  absorption  coefficient  at 
equilibrium  as: 

he  \ 

~kT(J) 

For  the  two-temperatures  model,  we  use  Eq.  89  and  the  Kirchhoff’s  law  at  Tve  to  estimate  the 
emission  coefficient. 


(89) 


=  Nmoi(T )S 


1  —  exp 


2.2.4  Photodetachment 

The  photodetachment  of  N_,  O-,  C_  which  is  described  by: 


A  +  hca  — >  A  +  e  (90) 

is  far  to  be  fully  understood.  The  available  theoretical  and  experimental  data  show  a  certain 
scatter.  According  to  the  literature,  the  major  contributions  correspond  to  the  photodetachment 
from  the  3P  term  of  the  ground  configuration  of  N_,  from  the  2P  term  for  O-  and  from  the  4S 
term  for  C-.  The  absorption  coefficient  has  been  expressed  as: 

he 
kT 

where  Nt,  the  population  of  the  above  identified  level  of  the  negative  ion,  is  estimated  through 
a  Saha  equation  from  neutral  atom  and  free  electron  populations  and  a  temperature  equal  to  T 
for  the  equilibrium  situation  and  to  Tve  for  the  two-temperatures  model.  The  cross  sections  Ej 
have  been  obtained  by  averaging  available  data  for  N_  and  O-  [10]  and  from  calculations  for 
C“  [73].  The  associated  emission  coefficients  are  calculated  using  the  Kirchhoff’s  law. 


=  [NiEi 


1  —  exp 
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2.3  Free-free  transitions 

The  free-free  emission  (Bremsstrahlung  emission)  is  observed  when  an  electron  is  slowed  down 
by  a  collision  with  an  ion,  atom  or  molecule.  The  energy  of  the  emitted  photon  is  equal  to 
the  loss  of  kinetic  energy.  The  cross  section  of  the  process  depends  on  the  relative  velocity  of 
colliders  which  can  be  taken  equal  to  the  electron  velocity.  As  the  electron  velocities  are  assumed 
to  follow  a  Maxwell-Boltzmann  distribution  at  Te,  the  bremsstrahlung  emission  coefficient  and 
the  inverse  bremsstrahlung  absorption  coefficient  are  linked  by  the  Kirchhoff’s  law  (Eq.  4)  with 
a  Planck’s  function  (Eq.  5)  at  Te. 


2.3.1  Electron-ion  interaction 


The  free-free  transitions  in  the  field  of  ions  are  taken  into  account  through  the  classical  Kramers 
formula  [74]  corrected  by  a  Gaunt  factor  g{a.  Te): 


■I,  =  9(<T,re)|  (gjjG 


2  a2e6  /  he 


mec° 


exp  (  cr  )  NeNi 


(92) 


where  Ne  and  Nt  are  the  densities  of  electrons  and  ions  respectively,  me  is  the  electron  mass,  a 
is  the  ion  charge  and  e  =  yf=- 

The  Gaunt  factors  of  Stallcop  and  Billman  [75]  have  been  used  above  11604  K,  those  of  Menzel 
and  Pekeris  [76]  below  11604  K. 


2.3.2  Electron-atom  interaction 


For  free- free  transitions  involving  atoms,  we  used  for  O  and  N  the  calculated  inverse  Bremsstrahlung 
cross  sections  Qa(a,Te)  of  Mjolsness  and  Ruppel  [77]  who  give  an  analytical  representation  of 
their  Hartee-Fock  calculation  for  the  temperature  range  4000-12000  K  and  the  wave  number 
range  4800-24000  cm-1,  and  for  C  the  ones  of  Geltman  [78]  which  cover  the  temperature  range 
500-20000  K  and  the  wavenumber  range  1000-20000  cm-1.  The  absorption  coefficient  is  given 
by: 


hb(j 


[Qa(a,Te)NeNat] 


1  —  exp 


(93) 


where  Ne  and  Nat  are  the  electron  and  atom  densities. 


2.3.3  Electron-diatom  interaction 

Free-free  transitions  have  been  included  for  N2  and  O2  using  the  cross  sections  of  Kivel  [79]  which 
cover  the  temperature  range  3000-15000  K  and  the  wavelength  range  0.3-4. 8  /mi.  These  cross 
sections  have  been  extrapolated  below  3000  K  and  below  2083  cnr1.  The  absorption  coefficient 
is  given  by  a  relation  similar  to  Eq.  93. 

The  total  absorption  coefficient  due  to  bound-free  and  free-free  transitions  for  an  atmospheric 
air  plasma  is  presented  on  Fig.  10  for  two  temperatures,  2000  K  and  15000  K.  The  contributions 
of  the  different  processes  are  presented.  The  preponderant  contributive  radiative  processes  are 
indicated. 


2.4  Validation 

An  effort  has  been  made  to  validate  the  selected  and  calculated  data  by  systematically  comparing 
the  results  with  available  measurements  or  predictions.  The  comparisons  have  been  made  when 
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Figure  10:  Continuum  absorption  for  an  atmospheric  air  plasma  at  2000  K  (left)  and  15000  K 
(right)  [10]. 


possible  for  a  given  process  (radiative  lifetimes  for  diatomic  electronic  systems  for  example)  and 
also  for  more  global  data  [9,  10,  11,  12], 

For  plasmas  at  local  thermodynamic  and  chemical  equilibrium,  energy  level  populations  are 
calculted  from  Boltzmann  distribution  and  the  chemical  composition  from  the  general  equilib¬ 
rium  laws  (Guldberg-Waage  equation  for  dissociation  and  Saha  equation  for  ionization),  and 
from  perfect  gas  law  and  electro-neutrality.  The  calculations  of  the  required  partition  functions 
at  LTE  or  for  a  multi-temperature  thermodynamic  state  are  described  in  [54]. 

Figure  11  presents  for  an  LTE  air  plasma  at  1  atm  the  optically  thin  source  strength  defined 


as: 


Sr 


KaI°{T)  da. 


(94) 


a i  and  02  correspond  the  experimental  cutoff  wave  numbers.  Our  calculated  source  strength 
corresponds  to  cutoff  wave  numbers  of  1000  cm-1  and  50000  cm”1  respectively.  These  values 
correspond  to  the  different  experimental  cutoffs,  except  the  ones  of  Nerem  and  Devoto  for  which 
(72  is  equal  to  58820  cm”1.  An  overall  good  agreement  is  observed  up  to  14000K. 


Spectral  comparison  has  been  made  with  the  emission  data  obtained  by  Laux  [80]  with  an 
ICP  torch  working  with  air  at  atmospheric  pressure.  Figure  12  presents  the  experimental  and 
calculated  absolute  emission  intensity  along  a  diameter  of  the  plasma  torch.  The  temperature  is 
about  7600  K  at  the  center.  It  must  be  noted  that  a  small  background  has  to  be  removed  from 
the  experimental  spectrum  [5].  Differences  in  the  UV  NO  bands  can  also  be  observed,  but  the 
overall  agreement  is  good,  especially  with  the  latest  measurements  given  in  Ref.  [5]. 


These  comparisons  between  experiments  and  predictions  in  the  equilibrium  regime  are  also 
relevant  because  the  problems  of  nonequilibrium  chemistry  and  non-Boltzmann  state-populations 
do  not  add  complexity  in  the  theoretical  modelling. 

Unfortunately  such  data  were  not  available  for  CO2  plasmas.  We  have  thus  developed  an 
emission  experimental  set-up  based  on  a  5  KW  microwave  torch  at  atmospheric  pressure  [86]. 
The  emission  spectra  of  CO2  and  CC>2(97%)-N2(3%)  have  been  recorded.  An  Abel  inversion  pro¬ 
cedure  allows  to  deduced  local  emissivities.  A  careful  analysis  of  the  emission  spectra  in  the  UV 
and  visible  ranges  has  been  performed  and  has  shown  that  the  plasma  is  at  local  thermodynamic 
equilibrium  for  molecules  with  a  slight  depletion  of  observed  atomic  levels.  The  predominant 
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Figure  11:  Optically  thin  emission  of  air  plasmas  at  atmospheric  pressure  in  the 
spectral  range  1000  -  50000  crn^1.  Comparison  between  our  calculations  and 
available  experimental  or  numerical  data,  a:  [80],  b:  [81],  c:  [82],  d:  [83],  e:  [84], 

f;  [85]- 


Figure  12:  Comparison  between  experimental  [80]  and  calculated  emission  intensity  along  a 
diameter  of  the  air  plasma  torch. 
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observed  systems  are  CN  violet,  CN  red,  and  C2  Swan.  Significant  contributions  are  also  found 
from  O2,  CO+  and  NO  systems.  Figure  13  presents  a  comparison  between  measured  and  cal¬ 
culated  emissivities  at  the  center  of  the  C02(97%)-N2(3%)  plasma.  The  molecular  temperature 
at  the  center  of  the  plasma  is  about  5500K.  For  the  pure  C02  plasma,  it  raises  to  about  6500K. 


400 
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600 


700 


800 


Figure  13:  Comparison  between  measured  and  calculated  absolute  emissivities  at  the  center  of 
the  atmospheric  C02(97%)-N3(3%)  plasma  (from  [86]). 


Figure  14  presents,  in  a  spectral  range  centered  on  the  CN  red  system,  the  experimental 
spectrum  and  the  best  adjustment.  A  systematic  study  of  the  sensitivity  to  the  temperature  has 
been  carried  out  by  calculating  the  root  mean  square  error.  Its  evolution  with  the  temperature 
is  given  in  the  insert  of  Fig.  14. 


2.5  Contributions  of  the  various  mechanisms  at  equilibrium 

We  analyse  in  this  section  the  relative  contributions  of  the  various  emission  mechanisms  to  total 
radiation  at  different  temperatures  for  air  and  CO2/N2  plasmas.  We  limit  our  discussion  to 
media  at  chemical  and  thermal  equilibrium  at  temperature  T  and  at  atmospheric  pressure.  The 
optically  thin  emission  strength  Sr  =  J KaI®(T)  da  is  shown  on  Figs.  15  and  16  for  air  plasmas 
(79%  N2,  21%  O2)  as  a  function  of  temperature.  Sr  is  limited  in  Fig.  15  to  the  spectral  range 
<7i=1000  cm-1  -  <72=50000  cm”1,  easily  accessible  experimentally,  while  the  emission  shown 
on  Fig.  16  includes  the  vacuum  UV  contribution  up  to  <72=150000  cm”1.  It  is  shown  that 
the  emission  is  dominated  by  diatomic  line  radiation  below  6000  K  and  by  atomic  lines  above 
8000  K.  Continuum  contribution  to  the  optically  thin  radiation  seems  to  be  relatively  small. 
However,  the  predominance  of  VUV  atomic  lines  at  high  temperature  is  generally  damped  in 
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Figure  14:  Comparison  between  measured  and  calculated  emissivities  fro  the  Av=3,  4  and  5 
bands  of  the  CN  red  system  at  the  center  of  the  atmospheric  C02(97%)-N2(3%)  plasma,  (from 
[86]). 


real  applications  due  to  the  strong  self- absorption  of  these  lines.  The  continuum  radiation  is 
then  no  longer  negligible.  To  outline  the  effects  of  absorption,  we  plot  on  figure  17  the  total 
plasma  emission,  as  corrected  from  self-absorption  over  a  characteristic  distance  R: 

POO 

Srs(R)  =  /  exp (-KaR)  da,  (95) 

Jo 

for  various  values  of  R.  At  15000  K  for  instance,  the  emission  is  reduced  by  one  order  of 
magnitude  for  R  as  small  as  1  mm. 


Figure  15:  Optically  thin  emission  of  air 
plasmas  at  equilibrium  at  atmospheric  pres¬ 
sure  between  1000  and  50000  cm-1. 


Figure  16:  Optically  thin  emission  of  air 
plasmas  at  equilibrium  at  atmospheric  pres¬ 
sure  between  1000  and  150000  cm-1. 
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Figure  17:  Volumetric  emission  of  air  plasmas  at  atmospheric  pressure  corrected  from  self¬ 
absorption  on  a  characteristic  distance  R. 


Temperature  (K) 


Figure  18:  Contributions  of  CO2  IR  bands,  diatomic  bound-bound,  atomic  bound-bound,  and 
continuum  radiation  to  the  optically  thin  emission  of  CO2  (97%)-N2  (3%)  plasma  at  atmospheric 
pressure  (adapted  from  Ref.  [12]). 
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Figure  19:  Details  of  the  relative  contributions  of  diatomic  and  atomic  bound-bound  radiation 
to  the  optically  thin  emission  of  CO2  (97%)-N2  (3%)  plasma  at  atmospheric  pressure  (adapted 
from  Ref.  [12]). 

For  C02(97%)-N2(3%)  mixtures  at  atmospheric  pressure,  Figure  18  shows  that  the  total 
optically  thin  emission  iis  first  dominated  by  CO2  infrared  radiation  (mostly  the  v 3  bands  near 
4.3  /mi)  below  4000  K,  then  by  diatomic  line  radiation  between  4500  and  6500  K,  and  finally 
by  atomic  line  radiation  for  higher  temperatures.  Here  also,  continuum  radiation  seems  to  be 
quite  weak,  but  VUV  atomic  lines  are  very  intense  and  are  strongly  self-absorbed  in  practical 
applications.  The  detailed  contributions  of  diatomic  and  atomic  bound-bound  transitions  are 
shown  separately  (without  CO2  radiation)  in  Fig.  19.  At  low  temperature,  CO  infrared  system 
is  the  most  contributing,  followed  by  CO  fourth  positive  system  between  4500  and  6500  K,  and 
then  by  carbon  atomic  lines.  It  should  be  noted  however  that  CO  IR  and  fourth  positive  systems 
are  likely  to  be  strongly  self-absorbed  in  practical  applications  as  discussed  in  Sec.  3.2. 


3  Application  to  earth  and  Mars  entries 

In  addition  to  the  spectroscopic  properties  described  in  the  first  part  of  this  paper,  the  com¬ 
putation  of  radiative  fluxes  in  entry  applications  requires  the  knowledge  of  the  entire  fields  of 
population  densities  of  the  emitting  states  of  each  species.  For  molecules  for  instance,  the  pop¬ 
ulations  of  all  electronic,  vibrational,  and  rotational  levels  are  required.  The  populations  are 
generally  calculated  by  neglecting  the  interaction  with  radiation  field  or  by  using  approximate 
escape  factors  in  some  collisional-radiative  models.  The  third  part  of  the  problem  concerns  the 
solution  of  radiative  transfer  problem  including  spatial,  directional  and  spectral  integrations. 

We  have  combined  our  spectroscopic  data  and  raditive  transfer  solvers  with  the  aero-thermo- 
chemical  fields  predicted  by  other  groups  to  compute  wall  radiative  fluxes  and  volumetric  radia¬ 
tive  powers  in  some  Earth  and  Mars  entry  applications. 
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3.1  Earth  entry 

The  most  widely  considered  data  for  radiation  model  validation  are  from  the  Firell  flight  exper¬ 
iment  [87].  The  database  has  been  used  to  calculate  the  radiative  intensities  at  the  stagnation 
point  for  4  points  of  Fire  II  re-entry  trajectory  [88].  The  aerothernral  fields  used  to  calculate  the 
different  spectra  are  those  of  Ref.  [89].  They  have  been  computed  using  the  chemical  kinetic 
model  of  Park  and  a  two-tenrperature  model:  Trot  for  the  translation  of  heavy  particles  and 
diatomic  molecule  rotation,  Tv  for  diatomic  vibration.  We  assume  that  electronic  levels  and 
free  electrons  are  in  equilibrium  at  Tv,  which  is  written  Tve  in  the  following.  Figures  20  and  21 
gives  the  temperatures  and  concentration  fields  along  the  stagnation  line  for  the  the  flight  time 

1642.66  s. 

Figure  22  shows  the  emission  and  absorption  coefficients  in  the  shock  layer  for  Firell  flight  time 


Figure  20:  Temperature  fields  along  the 
stagnation  line  for  the  trajectory  point 
1642.66  s  (from  [89]). 


Figure  21:  Concentration  fields  along  the 
stagnation  line  for  the  same  trajectory 
point  (from  [89]). 


1642.66  s  which  corresponds  to  the  maximum  radiative  heating.  The  most  important  mecha¬ 
nisms  correspond  to  atomic  lines  in  near  IR  and  VUV  regions,  N2  VUV  systems,  N^"  systems, 
and  photoionization.  The  contribution  of  VUV  represents  50%  of  the  total  incident  intensity  at 
the  wall.  It  can  also  be  seen  that  absorption  must  be  taken  into  account  in  the  visible  range  and 
especially  in  the  VUV  range.  Neglecting  absorption  in  the  calculation  leads  to  an  overestinration 
of  the  incident  flux  by  a  factor  of  55. 

Our  calculations  show  a  good  agreement  with  the  data  of  the  two  radiometers  (0.3-  0.6  fim 
and  0.2-4  /xnr)  for  the  trajectory  points  1642.66  s  and  1637.5  s.  The  calculations  overpredict 
the  radiative  intensities  for  the  trajectory  points  1634  s  and  1648  s.  Figure  23  compares  the 
measured  spectral  intensity  with  our  calculations  for  the  trajectory  point  1648  s.  A  few  molec¬ 
ular  systems  can  be  identified  (an  experimental  calibration  problem  explains  the  shift).  Inverse 
photodetachment  is  found  to  contribute  significantly  in  our  calculations.  This  suggests  to  reex¬ 
amine  the  modelling  of  this  process  to  decrease  the  associated  uncertainties. 

Finally,  figure  24  presents  the  cumulated  flux  incident  on  the  wall  for  the  trajectory  point 

1642.66  s.  It  points  out  the  necessity  to  take  into  account  the  effect  of  chemical  non-equilibrium 
on  O2  inverse  photodissociation 


3.2  Mars  entry 

As  there  are  no  direct  experimental  data  for  Mars  entries,  calculations  have  been  carried  out  in 
the  framework  of  testcase  TC3  of  the  ESA/CNES  radiation  working  group  to  simulate  radiative 
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Figure  22:  Spectral  emission  and  absorption  coefficients  for  a  point  in  the  shock  layer 
(Trot=16700  K  and  T^e=10250  K).  Flight  time  1642.66  s  (from  [70]). 


Wavelength  pm) 

Figure  23:  Spectral  incident  intensity  on 
the  stagnation  point  for  the  trajectory 
flight  time  1648  s  (from  [88]). 


Figure  24:  Cumulated  incident  intensity  on 
the  wall  stagnation  point  for  the  trajectory 
time  1642.66  s  (from  [70]). 
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fluxes  in  the  forebody  and  in  the  rear  region  of  an  axially  symmetric  vehicle  [90,  91]. 

Different  aerothermal  fields,  assuming  local  thermal  equilibrium  were  provided  for  this  test- 
case.  Figure  25  shows  for  instance  the  multidomain  temperature  field  around  the  vehicle  as¬ 
suming  a  pure  CO2  atmosphere.  We  carried  out  line  by  line  calculations  of  the  incident  flux 
on  different  points  of  the  vehicle.  Some  results  for  the  forebody  region  are  shown  in  Fig.  26 
and  the  computed  fluxes  in  the  rear  region  are  presented  later  as  an  application  of  statistical 
narrow-band  models. 

One  of  the  most  important  conclusion  of  these  calculations  was  that  the  incident  flux  is 
dominated  by  infrared  radiation  (mostly  from  CO2)  even  for  the  stagnation  line  where  the  tem¬ 
perature  exceeds  6000  K  in  the  shock  layer.  Figure  26  shows  indeed  that  the  intense  emission 
due  to  the  CO  fourth  positive  sytem  in  the  VUV  is  strongly  self  absorbed.  Its  dominant  contri¬ 
bution  to  the  optically  thin  emission  of  CO2  and  CO2/N2  plasmas  in  the  considered  temperature 
range  (see  e.g.  Figs.  18  and  19)  is  thus  strongly  attenuated  by  self-absorption.  This  result  also 
outlines  the  need  of  validation  of  CO2  spectroscopic  databases  above  3000  K. 


Temperature  (K) 


Figure  25:  Prescribed  temperature  field  around  the  simu¬ 
lated  vehicle  entering  Mars  atmosphere  [90]. 


4  Approximate  models 

The  use  of  a  line  by  line  calculations  to  perform  2D  or  3D  radiative  transfer  computations  is 
very  costly  since  it  requires  to  solve  the  radiative  transfer  equation  (eq.  1)  with  a  sufficiently 
fine  spectral  resolution  to  correctly  describe  the  transport  of  radiation  in  the  fine  structures  of 
the  spectra,  namely  in  lines.  This  motivates  the  development  of  approximate  models  which  do 
not  have  the  high  resolution  structures  but  provide  accurate  predictions.  Over  the  past  few 
decades,  such  models  have  been  developed  to  study  infrared  radiative  transfer  [92]  generally  for 
combustion  applications  or  atmospheric  propagation.  However  for  higher  temperatures,  such 
developments  are  scarce  [93]  and  to  our  knowledge  there  is  no  model  able  to  correctly  predict 
radiative  transfer  (emission  and  absorption)  in  a  non-equilibrium  plasma.  We  present  here  a 
model  specific  for  high  temperature-low  pressure  applications  at  thermal  equilibrium  which  has 
been  developed  to  estimate  the  flux  in  the  wake  of  a  an  axially  symmetric  body  entering  Martian 
atmosphere  [90]. 

4.1  Statistical  Narrow-Band  calculations  for  C02-CO  mixtures 

The  Statistical  Narrow  band  model  of  Mayer  and  Goody  [94]  has  been  used  to  model  the  infrared 
radiative  properties  of  CO-CO2  mixtures  between  500  and  5000  cm-1.  The  model  assumes 
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Figure  26:  Spectral  distribution  of  the  incident  radiative  flux  on  the  stagnation  point 
of  the  MSRO  simulation  (testcase  TC3  of  the  ESA/CNES  radiation  working  group 
[91]). 


random  line  positions.  The  transmissivity  of  a  uniform  column  of  length  £,  averaged  over  a 
spectral  band  of  width  A  a  is  expressed  as: 


ra  = 


A  G 


Laxl P(-^)d— p(-t) 


(96) 


where  5  =  A  a/N  is  the  mean  spacing  between  the  N  absorption  lines  located  in  A  a  and  W  is 
the  mean  black  equivalent  width  of  these  lines  defined  by  : 


_  i  _  r+o o 

w  =  jj  /  [l  -  exp  (~k1J)]  dm  (97) 

lines  i 


where  nla  is  the  absorption  coefficient  of  line  i. 

For  the  considered  application,  the  temperature  range  was  limited  to  4000  K  ,  and  the  pressure 
is  of  about  100  Pa.  The  line  profile  has  been  taken  of  Voigt  type  for  line  by  line  calculations 
but  assumed  in  the  model  of  Doppler  type  characterized  by  the  half-width  at  half-maximum 
7£>.  Among  the  different  line  intensity  distributions  which  have  been  tested  the  exponential 
distribution  has  been  selected.  It  yields: 


—  =  1Rf( 

5  5  V  7 d/S  )  ’ 


(98) 


where  u  is  the  pathlength  pa£,  pa  being  the  partial  pressure  of  the  absorbing  gas  and  E(y)  is 
defined  by: 
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(99) 


The  expression  98  introduces  two  parameters  S/S  and  7 d/5  which  have  been  determined  in 
narrow  bands,  25  cm-1  wide,  as  a  function  of  the  temperature  for  each  molecule.  S/S  corresponds 
to  the  mean  absorption  coefficient  divided  by  pa  in  the  weak  absorption  limit.  7 d/S  is  adjusted 
by  fitting  the  curve  of  growth  of  W/5  versus  u  obtained  from  line  by  line  calculations  for  optical 
paths  between  5  Pa.m  and  5000  Pa.m.  We  checked  that  the  mixture  transmissivity,  averaged 
over  A  a,  may  be  calculated,  with  an  excellent  accuracy,  as  the  product  of  CO  and  CO2  averaged 
transmissivities.  A  further  approximation  (Curtis-Godson  approximation  [92])  is  required  to 
compute  the  transmissivity  of  non-uniform  columns. 

The  incident  radiation  intensity  on  the  wall  along  a  ray  of  total  length  L,  averaged  overAu,  is 
given  by: 

W  =  Jo  /°(T(S))  JUr(S,L)ds  (100) 


4.2  Incident  radiative  flux  on  a  body  entering  Martian  atmosphere 

The  IR  incident  radiative  intensity  has  been  calculated  using  line  by  line  or  the  narrow-band 
model  for  the  aerothermal  fields  presented  in  section  3. 

A  ray  tracing  method  has  been  used  to  calculate  the  incident  flux  on  different  points  of  the 
vehicle  wall  with  900  discrete  solid  angles.  Due  to  the  low  pressure,  a  discretization  of  5  10~3 
cnr1  has  to  be  used  for  the  line  by  line  calculation  which  leads  to  solve  around  600  000  times  the 
RTE.  The  CPU  time  ratio  between  line  by  line  and  statistical  narrow-band  models  is  typically 
about  300.  The  good  agreement  between  the  two  simulations  is  illustrated  on  Fig.  28.  The 
third  curve  of  this  figure  corresponds  to  the  weak  absorption  limit  and  is  equivalent  to  the  use 
of  a  mean  absorption  coefficient  or  to  the  smeared  rotational  model.  By  neglecting  the  high 
resolution  spectral  correlations,  this  approximation  leads  to  an  overestinration  of  wall  fluxes  . 


Figure  27:  Calculation 
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Figure  28:  Comparaison  between  LBL  cal¬ 
culations,  full  SNB  results  and  SNB  results 
using  the  weak  absorption  limit  for  the  dif¬ 
ferent  basis  points  [90]. 


An  approximate  model,  able  to  deal  with  optically  thick  systems  of  diatomic  electronic  states 
under  nonequilibrium  conditions,  is  under  development.  In  the  same  manner,  such  model  for 
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CO2  under  vibrational  nonequilibrium  conditions  would  be  suitable  for  Mars  entries. 
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